Radiatively induced leptogenesis in a minimal seesaw model 
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We study the possibility that the baryon asymmetry of the universe is generated in a minimal 
seesaw scenario where two right-handed Majorana neutrinos with degenerate masses are added to 
the standard model particle content. In the usual framework of thermal leptogenesis, a nonzero 
CP asymmetry can be obtained through the mass splitting induced by the running of the heavy 
Majorana neutrino masses from their degeneracy scale down to the seesaw scale. Although, in 
the light of the present neutrino oscillation data, the produced baryon asymmetry turns out to be 
smaller than the experimental value, the present mechanism could be viable in simple extensions of 
the standard model. 



I. INTRODUCTION 

The outstanding advances witnessed by experimental 
cosmology have brought us to a new era where an un- 
precedented precision has been achieved in measuring 
several cosmological parameters. Of particular impor- 
tance is the value of the baryon-to-photon ratio, which is 
determined to be [1] 



Vb = 6.1 



+0.3 



10 



-10 



(1) 



from the latest measurements of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite. Several 
studies have been aimed at the explanation of this small 
but non-vanishing quantity. Among them, leptogene- 
sis [2] has become the most attractive mechanism to gen- 
erate the baryon asymmetry, not only because of its sim- 
plicity, but also due to its possible connection with low 
energy neutrino physics [3-6]. Indeed, it is conceivable 
that the physical effects responsible for the existence of 
suppressed neutrino masses may simultaneously play a 
crucial role in the thermal history of our universe. 

From the theoretical point of view, the addition of 
heavy neutrino singlets to the standard model particle 
content is an elegant and economical way to account for 
small neutrino masses via the seesaw mechanism. In this 
framework, the simplest possibility which can lead to a 
successful baryogenesis relies on the out-of-equilibrium 
L-violating decays of heavy Majorana neutrinos (thermal 
leptogenesis) . The excess of lepton number produced in 
these decays is then partially converted into a baryon 
asymmetry by the (B + L)-violating sphaleron interac- 
tions, which are in thermal equilibrium for temperatures 
10 2 GeV <T< 10 12 GeV. 

The heavy Majorana neutrino masses are rather un- 
constrained in the sense that they can range from a few 
TeV up to 10 16 GeV, depending on the model consid- 
ered. Moreover, their spectrum may be either hierarchi- 
cal, quasi-degenerate or even exactly degenerate. If the 
heavy Majorana neutrinos have a hierarchical mass spec- 
trum, a lower bound of 10 8 — 10 9 GeV can be inferred 
on the leptogenesis scale [7]. On the other hand, the 
resonant enhancement of the leptonic CP asymmetries 
through the mixing of two nearly degenerate heavy Ma- 



jorana neutrinos [8, 9] can lower their mass scale up to 
TeV energies [9, 10]. 

Due to the proliferation of free parameters in the high 
energy neutrino sector, the link between leptogenesis and 
low energy phenomenology can be only established in 
a model dependent way. In the framework of the SM 
extended with three heavy Majorana neutrinos one has 
eighteen parameters in the high-energy neutrino sector 
which have to be confronted with the nine low-energy pa- 
rameters of the light neutrino mass matrix. A minimal 
version of the seesaw mechanism, based on the existence 
of just two heavy neutrino states, reduces the number 
of parameters from eighteen to eleven. But even in this 
case, additional assumptions (such as the introduction of 
texture zeros in the Dirac Yukawa matrix) are required 
to completely determine the high energy neutrino sector 
from the low energy neutrino observables [11]. The effec- 
tive light neutrino mass matrix is in this case denned by 
seven independent quantities since one of the neutrinos 
is predicted to be massless. 

In this paper we assume an exact degeneracy of the 
heavy Majorana neutrinos at a scale which is higher than 
their decoupling scale. We show that, within the mini- 
mal seesaw scenario, the heavy neutrino mass splitting 
induced by renormalization group effects can lead to ap- 
propriate values of the CP asymmetries relevant for the 
computation of the cosmological baryon asymmetry in 
the thermal leptogenesis framework. This provides us 
with a simple explanation for the extremely small mass 
splitting which is typically required in this context. The 
compatibility of this framework with the presently avail- 
able low-energy neutrino data is also addressed. We 
show that the solar and atmospheric neutrino data puts 
strong constraints on this scenario and already excludes 
its minimal version. Since the radiatively induced lepton 
asymmetry turns out to be proportional to the charged- 
lepton r Yukawa coupling, an enhancement of this cou- 
pling could be sufficient to obtain the required baryon 
asymmetry. This can be achieved, for instance, in sim- 
ple extensions of the standard model with more than one 
Higgs doublet. In this case, the present framework could 
easily accommodate the low-energy neutrino data and 
lead to a successful leptogenesis. 
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II. SEESAW RECONSTRUCTION 



matrix 1 : 



In the framework of the SM extended with two heavy 
right-handed neutrino singlets, and before the decoupling 
of the heavy Majorana neutrino states, the lcptonic sec- 
tor is characterized by the Lagrangian 



-C = t L Yt£R<t> + D L Y v v R 



1)0 + -VrCMri/h + h.c. . 



(2) 



where ^l(r) and vl(r) are the left-handed (right- 
handed) charged-lepton and neutrino fields, respectively; 
0° is the neutral component of the SM Higgs doublet. 
The Dirac neutrino and charged-lepton Yukawa coupling 
matrices are denoted by Y v and Yg , while Mr stands 
for the 2x2 symmetric mass matrix of the heavy right- 
handed neutrinos. In the basis where Yt and Mr are 
diagonal, the neutrino sector is characterized, in general, 
by eleven parameters, namely the six moduli and three 
phases in the 3x2 Dirac neutrino Yukawa coupling Y v 
plus the two heavy neutrino masses, M\ and M 2 . This 
number may be further reduced by making assumptions 
about the structure of Y v [11-13]. 

Let us now assume that at a scale Ad the heavy neu- 
trinos are degenerate in mass, i.e. M\ = Mi = M, with 
M < A/5 . We remark that in the limit of exact heavy Ma- 
jorana neutrino mass degeneracy, CP is not necessarily 
conserved. Indeed, the non-vanishing of the weak-basis 
invariant [17, 18] 



Jx = M" 6 Tr \Y v Y?YtY}YfY2 ,Y?Y? 



1 3 



(3) 



which is not proportional to M| — Mf, would signal a 
violation of CP. 

It is convenient to define the 3x3 seesaw operator at 
Ar> , 



Y M _1 Y T — —Y Y T 
I v m R 1 v — j^*v x v ■ 



(4) 



An important feature of n is that, at any renormaliza- 
tion scale fi, one of its eigenvalues is zero leading to a 
massless neutrino at low energies. Therefore, this mini- 
mal seesaw scenario predicts a hierarchical neutrino mass 
spectrum and, as a consequence, a detection of neutrino 
masses with m, > yj Am, would automatically rule it 
out, independently of the seesaw parameters. 

For the cases we will be interested in, the decoupling 
scale of the heavy neutrino singlets is approximately M. 
Moreover, the RG effects in Y v and Mr can be neglected 
for the purposes of the seesaw reconstruction, although 
they will be crucial for leptogenesis. Thus, the seesaw 
operator at the scale M is well described by Eq. (4), 
apart from an overall factor of order one. 

In order to reconstruct the high energy neutrino sector 
in terms of the low energy parameters, we shall consider 
the following structure for the Dirac neutrino Yukawa 



Y v = yo 



Vi 

y3 2/4 
V5 ye 



(5) 



where J/1...6 are complex and y is chosen, without loss 
of generality, to be real. 

At low energies the effective neutrino mass matrix is 
given by 



A4 = m 3 U diag(0, p e ia , 1) , 



(6) 



where to 3 is the mass of the heaviest neutrino, a is a 
Majorana phase and p = mi/m^. The matrix U is 
the leptonic mixing matrix which can be conveniently 
parametrized in the form 



U 



c s c r 



s s c r 



SsCa C- s S a S r G C s C a SaSa$r& 



S a C r G 



iS 



iS 



(7) 



where Cj = cos 9j , Sj = sin 6j and S is the CP- violating 
Dirac phase. 

We will restrict our analysis to the case of neutrinos 
with normal hierarchy = mi < m 2 <C m 3 . In this case 



m 2 = ^Amf , m 3 ~ yj Ami , P - 



I Ami 
Ami ' 



(8) 



The case of inverted hierarchy can be analyzed in a sim- 
ilar way. In all our numerical estimates we shall use 
the best-fit values [15, 16] Am 2 s = 8.3 x 1CT 5 eV 2 and 
Ami = 2 - 2 x 1CT 3 cV 2 , so that p ~ 0.19. The atmo- 
spheric and solar mixing angles are taken as a = ir/A 
and tan 2 9 S = 0.37, respectively. The mixing angle 8 r is 
presently constrained to be | sin# r | < 0.2 at 95 % C.L.. 

The effective neutrino mass operator at the decoupling 
scale M can be found by running M. from the electroweak 
scale mz up to M . However, in the case when the spec- 
trum of the light neutrinos is hierarchical, these effects 
turn out to be irrelevant. In particular, the mixing an- 
gles and the ratio p are stable under the RG evolution. 
Therefore, one has 



m 3 [/diag(0,p e "M)£/ T ^— Y v Yj, (9) 

with v ~ 174 GeV. This relation leads to the following 
simple approximation for the reconstructed Dirac neu- 



1 For illustration we have chosen j/2 = 0. The same analysis can be 
extended to other one-texture-zero Dirac Yukawa matrices [13, 
14]. 
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trino matrix 
H = YjY„ ~ 

I x e 



Vl + 2m 2 cos a 



Q ia/2 



p + x 2 x [e*"'~ — pe 

-ia/2 _ pe ia/2j 1 + pa; 2 



ia/2] 



where 



y 



tan # r 



(10) 



(11) 



As expected, the assumption of one texture zero allows 
us to reconstruct the Dirac neutrino Yukawa matrix (5) 
in terms of low energy observables up to an overall factor 
which depends on the heavy Majorana neutrino mass M. 
To determine the latter, additional assumptions would be 
required. 



III. RADIATIVELY INDUCED CP 
ASYMMETRIES 

A non-zero leptonic asymmetry can be generated if and 
only if the CP-odd invariant 



J 2 = ImTr 



hm\m r m\e t m r 



= M X M 2 {M\ - M\) Im [H 2 2 



(12) 



does not vanish [9]. This requires not only the heavy 
neutrino mass degeneracy to be lifted, but also the non- 
vanishing of Im [H 2 2 ] at the leptogenesis scale M. Al- 
though the first condition is easily guaranteed by the 
running alone of M R from A^ to M, to achieve the sec- 
ond condition we must also include the quantum cor- 
rections due to the Dirac neutrino Yukawa matrix Y v . 
Indeed, the running of the right-handed neutrino mass 
matrix is governed by the renormalization group equa- 
tion (RGE) [19, 20] 



dM E 
dt 



= H 1 M R + M R H , t = 



In (m/Ad) 

16 7T 2 



(13) 



From this equation, and neglecting the running of Y v , one 
has in the leading- log approximation M R (t) oc 1 + (H T + 
H) t. It is then possible to show that the trace in Eq. (12) 
is a real quantity, which leads to J 2 — 0. Consequently, 
the CP-violating effects, to which the CP asymmetries 
are sensitive, vanish at the decoupling scale. 2 Thus, the 
corrections coming from the running of Y v from Ad to 
M must be taken into account. 

In the basis where M R is diagonal, and assuming that 
the charged-lepton Yukawa matrix Ye is diagonal, the 



This was first pointed out in Ref. [21]. 



evolution of the right-handed neutrino mass matrix and 
the Dirac neutrino Yukawa matrix is given at one-loop 3 
by [19, 20] 



^ = 2 Mi H u , 
dt ' 

dY v 
dt 



*-\&-\gl-\{Yrt-r*Yi) 



Y v + Y V R , 
(14) 



where T = 3Tr(Y u Yj) + 3Tv(Y d Yj) + Tr(Y £ Y/) + 
Tr(Y^Yjj'); Y u ^ are the up-quark and down-quark Yukawa 
matrices and gy,2 are the gauge couplings. The matrix 
R is antihermitian with 



Pll — P22 — 0, J?21 — — P12 



R 



2 + 5 



N 



12 



Re (#12) 



>N 



2 + 6 



N 



lm(H 12 ). (15) 



The parameter 



M 2 

5n = m[- 1, 



(16) 



quantifies the degree of degeneracy between M\ and M 2 . 

An important feature of the system (14)-(15) is the 
fact that, if the heavy Majorana neutrinos are exactly 
degenerate in mass (S N = 0) at a given scale A D , then 
the second RGE in (14) becomes singular, unless one 
imposes Ke(Hi 2 ) = 0. The latter condition can always 
be satisfied by performing an additional rotation in the 
Dirac Yukawa matrix, 



Yl = Y v O, O 



cos 9 sin 6 
— sin 9 cos 9 



(17) 



with the rotation angle 9 determined by the relation 

2Re(H 12 ) 



tan 29 = 

tl 22 — tin 

The rotated matrix H' = HO becomes 



H' 



Hn-A iIm(Hi 2 ) 
-ihn(H 12 ) H 22 + A 



(18) 



(19) 



where A = tan 6* Re (Hi 2 ). 

According to Eqs. (14), the running of the mass split- 
ting parameter Sn is determined by the equation 



N 



d5_ 
~dt 



= 2{l + 5 N ){H' 22 -H'n) 



(20) 



Given that Eqs. (10), (18) and (19) imply H' 22 - H' n ~ 
2/o(l — p), with yo defined in Eq. (11), the radiatively 



3 In general, the two-loop contribution to the running of Mjj enters 
in the calculation at the same level as the one-loop running of 
Y v . We have numerically verified that this contribution has a 
negligible effect on our results. 
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induced mass splitting at the decoupling scale M will be 
approximatively given by 



£L(l-p)ln(A D /M). 



(21) 



Moreover, even if at the degeneracy scale one has 
Re (H' 12 ) = 0, a nonvanishing real part will be generated 
by quantum corrections. The latter are easily estimated 
from Eqs. (14), and one finds 

Re(H[ 2 ) ~ M^ y 2 Re(eW 2^* 3)ln(AD/M) 7 



(22) 



where y T is the r Yukawa coupling and U is the neutrino 
mixing matrix defined in Eq. (7). 

The crucial ingredient for the computation of the 
baryon asymmetry in a leptogenesis scenario are the CP 
asymmetries generated by the interference of the one- 
loop and tree-level heavy Majorana decay diagrams. In 
the case when the heavy Majorana neutrinos are nearly 
degenerate in mass, these asymmetries are approximately 
given by [10] 




FIG. 1: Comparison of the analytical and numerical result for 
the CP asymmetries ei and £2 as functions of s r = sin 6 r . The 
analytical curves are defined in Eqs. (25), while the numerical 
ones were obtained by integrating the full set of RGEs from 
the heavy neutrino degeneracy scale Ad to the decoupling 
scale M. 



Im[^ 2 ] 
16 kSnH'jj 



4M 2 S 2 N 



where 



HIM 
8tt 



.7 = 1,2, (23) 



(24) 



are the tree- level decay widths 4 . We notice that the ex- 
pressions (23) for the CP asymmetries £j exhibit the ex- 
pected enhancement due to the mixing of two nearly de- 
generate heavy Majorana neutrinos [8, 9]. 

In terms of the low-energy neutrino parameters, we 
obtain 



£1 



3y 2 % /p(l + p)x 



647r(l-p)(p + 
x sin(a/2) [c s cos(d - 



: 2 - A) 

a/2) + ^/p~s 2 s xcos(a/2)] 



can be readily seen from Eqs. (25). This corresponds to 
s r = ps s ~ 0.1 and 



3y 2 c s (l+ P ) 

1287T (1 - p) 



-10- 



(26) 



The accuracy of these approximate expressions is 
shown in Fig. 1, where the CP asymmetries are plotted 
as functions of s r taking A D = 10 16 GeV, M = 1 TcV, 
5 = it/2, a = 7T and assuming y T = 0.01 in the analyti- 
cal estimates. The solid lines correspond to the full nu- 
merical integration of the RGE system, while the dashed 
ones refer to the approximate values given in Eqs. (25). 
The comparison of the curves shows that, for values of 
s r ^ 0.1, £2 is suppressed with respect to £1, in accor- 
dance with Eqs. (25). 



£2 _ p + x 2 - A 
£1 ~ 1 + px 2 + A ' 



(25) 



A = -(1-P) 



-1 + x 2 + \/l + 2x 2 cosa + x 4 



Taking for instance a ~ n and S ~ 7r/2, the CP asym- 
metry £1 reaches its maximum value for x — ^fp , as 



4 One may wonder whether finite temperature effects could 
modify the perturbative result. Since the thermally induced 
mass splitting is approximately given by M2(T) — M\(T) < 
j^Rc(Hl) 2 iT 2 /M [10], at T ~ M this contribution is very sup- 
pressed with respect to the zero-temperature mass splitting 5jv- 



IV. COSMOLOGICAL BARYON ASYMMETRY 

The out-of-equilibrium Majorana decays are controlled 
by the parameters 



Ki 



H(T = Mi) ' 



(27) 



where H(T) = 1.66gl /2 T 2 /M P is the Hubble parame- 
ter, g* ~ 107 is the number of relativistic degrees of 
freedom and Mp = 1.2 x 10 19 GeV is the Planck mass. 
Assuming that the entropy remains constant while the 
universe cools down from T ~ M to the recombina- 
tion epoch, one can estimate the baryon-to-photon ratio 
as rjB — — 10~ 2 (di E\ + ^2^2) , where di < 1 are dilu- 
tion factors which account for the washout effects. Since 



■5 




s 

r 



FIG. 2: Numerical result for the baryon-to-photon ratio t\b 
as a function of s r — sin^r, obtained by integrating the full 
set of RGEs from the heavy neutrino degeneracy scale Ad to 
the decoupling scale M. 



from Eqs. (10), (25) and (27) it follows that \e 2 /K 2 \ < 
\ei/Ki\, one expects 

~ -10- 2 d l£l . (28) 

We also note that K\ is independent of M and approx- 
imately given by 

44 (P + ^ 2 -A) (29) 
V 1 + 2x 2 cos a + x A 

In this washout regime, a simple estimate of the di- 
lution factor d\ can be obtained from the fit d\ ~ 
0.6 [ln( J ft'i/2)]- - 6 / J ft:i [22]. For a ~ tt, (5 ~ tt/2, the max- 
imal v alue of the baryon asymmetry is then attained for 
x ~ v / 3p/(l + p)/3 ^ 0.23 and s r ~ 0.05. From Eq. (25) 
we find e\ ~ — 8 x 10~ 7 . Moreover, since Eq. (29) im- 
plies K\ ~ 11, then rfi ~ 4 x 10~ 2 , using the above fit. 
Inserting these values into Eq. (28) one finally gets 

77^ ax ~ 3 x 10~ 10 , (30) 

which is by a factor two smaller than the observed baryon 
asymmetry [cf. Eq. (1)]. It is worth noticing that this 
result is weakly dependent (through the renormalization 
effects on j/ 2 ) on the heavy Majorana neutrino mass scale 
M, as can be seen from the approximate expressions 
given in Eqs. (25). In Fig. 2 we present the numerical 
estimate for the baryon-to-photon ratio t\b as a function 
of s r . 

We remark that a more accurate computation of the di- 
lution factors requires the solution of the full set of Boltz- 
mann equations. This question has been readdressed in 
recent works, where the effects on leptogenesis due to the 
AL = 1 processes involving gauge bosons [10, 23] and 
thermal corrections at high temperature [23] have been 



discussed. However, since the present neutrino data fixes 
the decay parameters Ki to be Ki ^> 1, the above effects 
turn out to be negligible and do not alter significantly 
the final value of the baryon asymmetry. 

V. CONCLUSIONS 

We have considered a minimal seesaw scenario with 
two heavy Majorana neutrinos which are degenerate in 
mass at a scale higher than their decoupling scale. We 
have shown that the heavy neutrino mass splitting, in- 
duced by radiative corrections, leads to non-zero CP 
asymmetries, and subsequently, to thermal leptogenesis 
even for heavy Majorana neutrino masses as low as a 
few TeV. However, we have seen that in the minimal ex- 
tension of the SM with the addition of only two heavy 
Majorana neutrinos, it seems difficult to reconcile the 
present solar and atmospheric neutrino data with the 
observed cosmological baryon asymmetry On the other 
hand, the analytical expressions (25) already suggest that 
these constraints can be easily relaxed, for e.g., in exten- 
sions of the SM with more than one Higgs doublet, since 
in this case the charged-lepton Yukawa couplings can be 
substantially enhanced. 

We also remark that the results obtained here depend 
on the specific structure of the Dirac neutrino Yukawa 
matrix. In particular, it would be interesting to disen- 
tangle the texture considered here from other seesaw tex- 
tures which can lead to the same low-energy observables. 
A more general analysis regarding this issue will be pre- 
sented elsewhere [14]. It would also be desirable to find 
a symmetry which naturally leads to degenerate heavy 
Majorana neutrinos at high scales and, simultaneously, 
to the required Dirac neutrino Yukawa couplings. 

Although low-energy (TeV-scale) leptogenesis is viable 
in the present framework, it requires tiny Dirac neutrino 
Yukawa couplings of the order of 10~ 6 — 10~ 7 , which 
makes it difficult to be testable at future colliders. On 
the other hand, the fact that the final baryon asymmetry 
weakly depends on the heavy Majorana neutrino scale 
allows one to lower this scale below the typical values 
(<~ 10 8 ~ 9 GeV) required in the standard thermal lepto- 
genesis scenarios. Thus, the potential problems related 
with the overproduction of relic abundances, which could 
jeopardize the successful nucleosynthesis predictions, are 
avoided. 

Note added: While this work was being revised, 
Ref. [24] appeared, where it is shown that in the super- 
symmetric minimal seesaw model the present mechanism 
is indeed viable and thermal leptogenesis is successful. 
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